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Abstract 



We compute the exact spectral density of random matrices in the ground state 
of the quantum hamiltonian corresponding to the matrix model whose double scaling 
limit describes pure gravity in 2D. We show that the non-perturbative effects are very 
large and in certain cases dominate the semi-classical WKB contribution studied in 
the earlier literature. The physical observables in this model are the loop averages 
with respect to the spectral density. We compute their exact ground-state expectation 
values and show that they differ significantly from the values obtained in the WKB 
approximation. Unlike the alternative regularizations of the nonperturbative 2D quan- 
tum gravity, based on analytic continuation of the Painleve transcendent, our solution 
shows no pathologies. 
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The discovery [|TJ of the double scaling limit of the matrix models of 2D Gravity and 
its relation || to the KdV hierarchy opened the way to the study of the nonperturbative 
phenomena in Gravity and String theories. These phenomena are supposed to be described 
by the so called string equation 

x = $>(2J + l)iMC] (1) 

where x is the scaling combination of cosmological constant and parameter 1/N of the genus 
expansion, ti are the mass parameters of the theory and Ri [C] = C l + derivatives are the 
Gelfand-Dikii differential functionals corresponding to the specific heat C(x) as the potential 
of the Schrodinger operator ( see e.g. |§ for details). 

The string equation is in general a differential equation for C(x) of the order correspond- 
ing to the highest / in the sum. As observed in Q, for odd / there is the natural boundary 
condition 



C(x) -> x 1/l x-^±oo (2) 

which uniquely determines the solution. 

For the even / case situation is more complicated. The above boundary condition would 
lead to a complex solution, corresponding analytic continuation of the divergent matrix 
integral ||. The attempt to define the even I case as the limit ti + \ — ► of the odd case || 
does not work, since the limit does not exist. 

One way to interpret this paradox is to claim the even case and in particular the I = 2 case 
of pure 2D Gravity inconsistent. That would be an exciting possibility : the nonperturbative 
creation of matter in Gravity. However, the matter here is not quite physical, as it does not 
correspond to the unitary multiplets of conformal field theory. The partition functions of 
the I > 2 models with higher genus are not positive 0. 

So, in fact, this interpretation implies that no 2D Gravity could exist at nonperturba- 
tive level. This is too serious a statement to make without investigating the alternative 
interpretations. 

One such interpretation, suggested by Marinari and Parisi along the lines of the gen- 
eral method by Greensite and Halpern ||, seems particularly appealing. They do not add 
any ad hoc terms to stabilize the equation, but rather modify the basic definition of the 
matrix model to make it meaningful for arbitrary potential. In Ref. || the spectral density 
of the resulting quantum hamiltonian was studied in the semi-classical WKB approximation. 
It turns out to be possible to go beyond the WKB approximation and to obtain the exact 
spectral density of this hamiltonian. The necessary framework was set up in Ref. [To|, in 



which the set of equations generalizing the string equation was derived and studied numer- 
ically. Quantum effects, invisible in the WKB approximation, turned out to be extremely 
important. 



The prescription of Refs. ||,|§ has attracted a considerable amount of attention, |TT 
In particular, there have been some interesting attempts to study the properties of 
the model numerically at finite N. The problem encountered in these studies is the same 
as with all numerical studies of the matrix models, namely the double scaling limit requires 
very large N, because of the N 1 ^ 5 dependence of physical observables. 



2 



The goal of the present paper is to apply the analytic and numerical methods of the 
note [|l(J to derive the exact spectral density of random matrices in the ground state of 
2D Gravity. As a result, we are able to demonstrate that the non-perturbative effects are 
very large and in certain cases completely dominate the semi-classical WKB contribution. 
The physical observables in this model are the loop averages with respect to the spectral 
density. We compute their exact ground-state expectation values and show that they differ 
significantly from the values obtained in the WKB approximation. 

We begin by discussing the method of Greensite and Halpern for the simplest example of 
a one dimensional integral. Consider a positive action S(x) depending on a single variable 
x: 

dxexp [— S(x)} (3) 



By definition, the average of an operator Q(x) in the action S is given by 



dxQ(x) exp[— S(x)] 

(Q)z = — r (4) 

dxexp [—S(x)] 



If we define 

ij = exp [-S(x)/2] /VZ (5) 
the expectation value can be written 

(Q)z = (HQ\A) (6) 



i/jq is always positive and therefore it is always possible |22 ,||,[f| to construct a quantum- 



mechanical hamiltonian H = —d 2 /dx 2 + V(x), such that i[)q is its ground-state wave-function 
with the ground-state energy equal to zero: 

H^ = 

-^(x) + V(x)Mx) = (7) 
V(x) = ip'o(x)/ijjo{x) 

It is instructive to rewrite this hamiltonian in a manifestly positive form, as a product of an 
operator and its hermitean conjugate, 

"00 is annihilated by the second factor and hence Hifj = 0. When the action S is unbounded 
from below, the average (f|) is ill-defined. Formally, (||) is still annihilated by (||), but ipo is 
now not normalizable and therefore cannot be the ground state of H. On the other hand, 
H is a positive hamiltonian by construction and therefore a normalizable ground state \I/o 
must exist. Clearly the true ground state must be different from ip and must have an 
eigenvalue e > 0. Consider the action S(x) = x 2 — g 2 x 4 . The average (||) is ill defined, but 
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the power-series expansion in g 2 exists and can be computed explicitly. Each order in the 
expansion in powers of g 2 involves only moments of x with a gaussian measure, 

00 (g 2 x 4 ) n 
exp[— S(x)] = exp(— x 2 ) ^ 



J dxQ(x)exp(-x 2 ) £ (g 2 x A ) n /n\ 

n=0 

/ dxexp(-x 2 ) (g 2 x 4 ) n /n\ 



(Q) l z = (9) 



n=0 

Any sensible definition of average with the action S(x) ought to reproduce the perturbation 
expansion (||) and it should reduce to @ for a positive action. The quantum-mechanical 
expectation value in the true ground-state \l/o of H satisfies both constraints, we therefore 
take it as the definition of the average with the action S: 

(Q) z = <*o|Q|*o> (10) 
This prescription can be applied to matrix models, 

e F = Z = J d(j)exp [-(3tvU k ((j))] (11) 

where is a hermitean N x N matrix and the critical potentials are given by 

- 1 dt 



U k {4>) 



t 



(l-t(l-t)0 2 ) 



(12) 



At large <fi the 4> 2k term dominates, with the coefficient (— l) fc+1 Jq 1 dtt k ^ 1 (l — t) k . Therefore 
the fc-even critical potentials are unbounded from below and it is necessary to provide a 
prescription for defining the average (0) beyond perturbation expansion. Following the 
work of Greensite and Halpern || and Marinari and Parisi |§, we adopt the prescription 
( [lOD . After all, the only existing justification of the matrix models as theories of gravity 
is perturbative, in the sense of the genus expansion. The quantum mechanical definition 
is a priori as good as the statistical one, but has the advantage of being guaranteed to 
make sense for all models. The extra bonus of this prescription is the relation to the one 
dimensional supersymmetric string theory ||, with dynamically broken supersymmetry in 
the even I case. For practical purposes the supersymmetry seems to be useless so far, but 
one may hope to relate it to the physical supersymmetric string theories. 

The simplest case of an unbounded potential in a matrix model is the cubic potential 

U (<t>) = Y~ 9 3 (13) 

In order to proceed, it is useful to change variables from the hermitean N x N matrices to 
their eigenvalues, {Ai, . . . , Ajv}: 

e F = Z = J\{d\ n(A, - \j) 2 exp [-/3tr U(X) ] (14) 

i i<j 
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The factor TT(Aj — A^) 2 comes from the Jacobian. In the following it is useful to include this 

i<j 

factor in the wavefunction. The integration measure is then simply / J^dAj. It follows that 

J i 

the effective wave function corresponding to ( |l4l) is 

^> eff =\{{X l -X J )^{X)/^Z (15) 

i<j 

The hamiltonian H is immediately obtained from (|3|) and (|TH): 
H = J2H t ; H i = ~ + V p g{X i ); V eff (X) = N(3 (gX - I) + ^ (A - gX 2 f (16) 



OX ! " " ' ' V 2 J 1 

(15) and ([TB]) describe the ideal Fermi gas of N non-interacting particles. The fermionic 
property is generic - it is due to the Jacobian of the transformation from the matrices <p to 
their eigenvalues. On the other hand, the fact that the effective fermions are free is only true 
when the original potential U{<p) is cubic in 0. 

The effective potential depends on the free variable A and on three parameters of the 
original matrix-model action: N, (3 and g, V e jj = V e jj{N, j3, g; A). We are interested in the 
critical properties of the theory in the double scaling limit, 

x = ^^; iV-oo; x~N° (17) 

with the critical point at x — 0. For small positive x the effective potential V e jj(N, [3, g; A) 
has a double-well shape. In the double scaling limit ( |17D the depth of the left well scales like 
~ iV 2 , while the depth of the right well scales like ~ x 3 ^ 2 ~ A^°. Despite this, it is the right, 
tiny well, which is responsible for all the interesting critical phenomena in the double scaling 
limit. All the cubic potentials of the type ( |T3D are in the same universality class. At the 
critical point the depth of the second well goes to zero, and the potential has an inflection 
point, instead of a double minimum. 

Since we are interested in having the critical point at x = 0, it is necessary to chose a 
value of g such that for x = 0, i.e. for (3 = N, the extremum is an inflection point: 



dV eff 



leading to 



= 0; 

A=A 



PVeff 



dX 2 



= (U 

A=A 



'12^3 ^ ^ 

A = (73 + l) ^3 (20) 

We are interested in the scaling properties of the Fermi energy, or ground state, of the 
hamiltonian fllCf ). The scaling properties are determined by the tiny right well, but the bulk 
of energy levels is in the left, bigger well. The Fermi energy, ep can be written as 

e F = e% + e s F {x) (21) 
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where the bulk part, e° F , is determined by the left well and e s F (x) is the ^-dependent scaling 
part that we are interested in. The depth of the left well ~ N 2 , while the depth of the right 
well ~ N°. In terms of magnitude, cf is completely dominated by e F . The hamiltonian 
(|16|) can only be solved by some approximation procedure, analytical or numerical. It is 
mandatory to first isolate and "magnify" the scaling part, otherwise any approximate result 
for e F will be dominated by e F and the x dependence of the fine-structure will be lost. 

To isolate the scaling part, we expand the effective potential in x and in A — Ao around 
the inflection point and obtain the scaling potential v (y): 

V eff (N, (3, g- A) = V eff (N=P, g; A=A ) + (3 A 'Wv{y) (22) 
v(y) = j-ey (23) 

where x is given by (|17]) and 

y = (A - \ )a^P^, e = xa~V\ a = ^2 (24) 

The critical properties of the theory are thus determined by the scaling hamiltonian 

h = -£~ 2 +v(y) (25) 

The cubic potential v(y) (not to be confused with the original cubic potential fll3|) ) is for- 
mally unbounded from below. This is not a problem, however, since v(y) is only meaningful 
in the scaling region and the full potential V e fj{\) is bounded. The strategy for computing 
the fine structure of the Fermi energy is then as follows. An exact numerical solution for the 
spectral density of the scaling hamiltonian (|25| ) can be obtained through the powerful meth- 
ods of Gelfand and Dikii [Q. The exact solution can then be compared with the spectral 
density obtained from the WKB approximation. The two differ in the scaling region only, 
since WKB is completely adequate in the large left well. Their difference converges fast, 
yielding the contribution of the scaling region. 

In order to obtain the density of states, we first write down the general equations of the 
Fermi-gas theory in one dimension. It is convenient to scale out (3 from the Hamiltonian, 
and introduce the resolvent, 

G(e,y) = (y\(h-e)- 1 \y} (26) 
The particle density p(e,y) is related to the imaginary part of the resolvent, 

p{e,y) = -ImG(e + iO,y) (27) 

7T 

and the spectral density v{e) is given by the integral of p 

POO 

v{e)= dyp{e,y). (28) 
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The normalization of v[e) is fixed by the fact that ([15]) and ([16]) describe N non-interacting 
fermions. We rescale the spectral density by TV and the equation for the Fermi energy e B is 
therefore 

1 = / deu(e). (29) 



<x> 



The WKB particle density is given by 



PwKB(e,y) = y (30) 

27i^e-v(y) 

with the corresponding WKB spectral density v WKB 

ryi 1 

vwKB{e) = dy (31) 

J -°° 2nJe — v(y) 

where y\ is the first root of e — v(y). v WKB {e) can be expressed in terms of elliptic integrals. 
When there are 3 real roots y\ < y 2 < y-s of v(y) — e = 0, from eqs. (3.131.1) and (8.112.1) 
of Ref. p5fl we obtain 



vwKB(e) = — - 1 K(p); p = J— — (32) 

^yy-i - yi V 2/3 - yi 

When there is one real root y\ and two complex-conjugate roots y\ = y%, one can use 



,WKB 



eq. (8.126.1) of Ref. |25j to transform (|32D into the form 

VWKB{e) = K(sin(0/2)); y 2 - Vl = rje^ (33) 

The "Fermi energy" corresponding to the WKB solution, e F /KB is defined by 

dev WKB {e). (34) 

Outside the scaling region the exact spectral density is equal to the WKB density, 

v(e) = v WKB {e) for e < ei ow , (35) 

where e\ ow is some large negative value of e. The WKB solution for the Fermi energy e F 7KB 
for e > was found in . It exactly coincides with the bottom of the second well, which in 
our normalization is 

e W F KB = -\eV\ (36) 



Combining (|29|), (|34"P, (|35|) and (0), we obtain the final implicit equation for the Fermi 
energy e F as determined by the scaling region: 



r e F KB r F 

/ deu WKB (e) = / deu(e). (37) 
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The exact solution for spectral density can be obtained from the Gelfand-Dikii equation [24 
for the resolvent G(e,y) 

- 2Gd 2 G/dy 2 + {dG/dyf + 4(u - e)G 2 = 1 (38) 

Under usual circumstances, when the parameters are not fine-tuned to magnify the scaling 
region, the G-D equation is rather useless, as the direct solution of the Schrodinger equation 
is simpler. However, in our case it is just what we need. Differentiating ( |3~8D with respect 
to y, dividing by G and taking the imaginary part, we obtain an ordinary third-order linear 
differential equation for the continuous particle density p(e, y) in the double scaling limit, 

p" = 2v'p + 4(u - e)p (39) 

At large \y\, the general asymptotic form of the solution for p, correct up to terms ~ 
can be written in terms of three integration constants cf, cf and cf: 



exp ^2 J ^dy^j + cf exp ^—2 J y^vdy) + cf 



y -> ±oo (40) 

To chose the proper boundary conditions, physical intuition about the system must be used. 
At y — > oo the potential v(y) grows large, so only the decaying exponential is left, while 
cf = cf = 0. This leaves one free parameter, cf , the overall normalization of p, which 
can be determined as follows. We start from the asymptotic solution ( f40|) at some large 



positive y = yo, with cf = cf = and some arbitrary initial value of cf. We then solve 
the differential equation fl39|) numerically, evolving down to large negative values of y, where 
v(y) < and y/v is complex. The solution there is of the form (0), with cj~, C2 , C3 7^ 0, i.e. 
it contains two oscillating exponentials, plus a powerlike term. The oscillations represent 



a pure quantum effect, invisible in WKB expansion. If we average py (e — v) over these 

oscillations, the c{ and c^ terms disappear and we should obtain according to the 
WKB solution. This means that the solution ought to be multiplied by a constant such that 
C3 = 1/27T. This is the missing normalization condition for the density. 

In practice averaging over the oscillations is rather tricky, as it involves delicate cancel- 
lations between the positive and the negative contributions. There is, however, a better way 
of extracting C3 from p. If p is given by (flCfr ), then, up to terms ~ 0(1/ \v\), 




c,=^^\p--^—\ (41) 



The prescription (|41f) has the advantage that it is local in y and only requires the knowledge 
of p", which is readily available in any code used for solving differential equations. The 
reliability of fl4lf) can easily be tested, by verifying that the result for C3~ is independent of y. 

We employed the differential equation solving routine ODE, described in detail in Ref . |2(J . 
The routine is very stable and extremely easy to use. It is based on a variable-step, variable- 
order Adams method (explicit linear multistep method). The variable-step feature is essen- 
tial, because of the rapid crossover from smooth to oscillatory behavior. 
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Once p{e,y) is known, v{e) is in principle given by (|28|). In practice, it is more efficient 
to solve for the integral of p(e,y): 



£(e,y) = -/ dyp(e } y); p(e, y) = <9£(e, y)/dy (42) 

J y 

which satisfies the 4-th order differential equation 

= 2v'i' + A(v - e)£" (43) 
yielding both p(e, y) and z/(e) at the same time, 

v{e) = - lim £(e,y) (44) 

y -» -oo 

The boundary conditions for the derivatives of £ are those for p and its derivatives. In 
addition, £(yo,e) = 0, where yo is large and positive. 

While solving eq. (fffi) for the various values of e, we store the values of p(e, y) on a two- 
dimensional grid in the [e, y] plane. These are to be used later in obtaining the distribution 
of the eigenvalues of a random matrix (see eq. ( fi^D below). 

The tail of the integral in (fE|) converges rather slowly, since p(e,y) ~ 1 2/ 1 3 / 2 f or large 
negative y. In that region, however, the integral of p(e, y) is very well approximated by the 
integral of Pwkb{^,v), eq. flSDP. The corresponding WKB tail r(y) can expressed in terms 
of elliptic integrals, which we calculated numerically using standard methods. 

The solution of the G-D equation for p at e = 1 and e = — |, which corresponds to the 
bottom of the second well is presented in Fig. 1. The quantum effects are huge! 

Once v[e) is known, we obtain the Fermi energy ep(e) from eq. fl37j). The result is shown 
in Fig. 2. 

We can now obtain the distribution £l e (y) of the eigenvalues of a random matrix in the 
double scaling limit. This is done by integrating p(e,y) with respect to e, up to e^(e): 

n e (y)= [ eF(t) dep(e } y) (45) 



In practice the integrand is obtained by interpolating p(e, y) from the previously stored two- 
dimensional grid in the [e,y] plane. The resulting Q e {y) is shown in Fig. 3. and compared 
with the corresponding WKB density, 



a WKB 



VY KB (y) = r"~ {e) dep WKB (e, y ) = ^e™ KB -v{y) (46) 

J —00 7T 

The contribution of the non-perturbative effects in physical observables can most easily 
be seen by computing loop averages for positive I: 

1 r°° 

W e (l) = — d V Q e ( V )e lv (47) 

Z/ e J — 00 

where _ 

drjQ 6 {ri) (48) 
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The integral in ( (4"T| ) converges fast on both ends. For large positive y, one is in the classically 
forbidden region, where Q e decreases faster than exponentially (see Fig. 3), and for large 
negative y the exponential factor in the integrand ensures exponential convergence. Thus 
in practice it turns out to be sufficient to take a finite upper limit of the integration at 
y = yright, i-e. the value of y at which fl e (y) becomes sufficiently small. 

Comparison of W e (l) with its WKB analogue, defined by eqs. d^),© with Q^ KB {y) 
instead of Q t (y), shows that the non-perturbative part of spectral density takes over the 
WKB part for large I. The results for W e (l) for < / < 6 are shown in Fig. 4. 

Note, that the exact loop average, unlike the WKB one, reaches the minimum, and then 
grows! This is a direct consequence of the fact that the scaling eigenvalue density fl e {y) is 
different from zero in the classically forbidden region y > 0. 

In physical terms, this striking phenomenon is explained as follows. The larger the 
boundary I of our two dimensional space, the more area it could encircle, and hence, the 
more handles could be attached to it. The number of surfaces with unlimited genus grows 
as exponential of the fifth power of the area 0, and a typical area grows as the square of 
the length of its boundary. So, we might expect a very fast growth of the loop average due 
to nonperturbative contributions. 

Another comment: in the literature, doubts were raised |27 , whether the Marinari-Parisi 
prescription (in the WKB approximation) could preserve the positivity of the loop average. 
We did not observe such pathologies. Moreover, we believe that with correct nonperturbative 
definition of the loop average the positivity is guaranteed. 

The point is that the spectrum of the random matrix in the Marinari-Parisi model is 
infinite, it covers the whole real axis. Our spectral integral for the loop average is manifestly 
positive definite. 

However, in the WKB approximation, studied in J271, there are several disconnected 
regions of the classically allowed motion with gaps in between. Our prescription is to add 
all these contributions with positive sign. Being rewritten as a contour integral, this would 
produce a set of loops encircling each allowed region anticklockwise, or a single contour going 
in the imaginary direction to the right of all the regions. 

For positive I one can close the contour in the left half plane, yielding the original spectral 
integral. For negative I one could close the contour in the right half plane, yielding zero, as 
there are no singularities of the integrand. 
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Figure Captions 

• Figure 1. The particle density p(e,y), solution of eq. (|39|), for e = 1 and e = — |. 
Dash-dotted line denotes the WKB solution. 

• Figure 2. The Fermi energy, eir(e), obtained from eq. ([37]). The diamonds denote the 
actual values computed, the continuous curve is plotted to guide the eye. Dash-dotted 
line denotes the Fermi energy in the WKB approximation, ep KB (e). 

• Figure 3. The scaling eigenvalue density £l e {y). The WKB approximation is plotted 

dash-dotted line. 

a) Q e {y) and n? KB (y), e = -1, 0, 1. 

b) AQ e (y) = ^ KB {y) - Q e (y), e = -1, 0, 1. 

• Figure 4. Loop average W e (l), eq. (|47|), for e = —1, 0, 1. The WKB approximation is 
plotted dash-dotted line. 
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